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A  meso-scale  unit-cell  based  material  model  for  a  prototypical  plain-woven  single-ply  flexible  armor  is 
developed  and  implemented  in  a  material  user  subroutine  for  use  in  commercial  explicit  finite  element 
programs.  The  main  intent  of  the  model  is  to  attain  computational  efficiency  when  calculating  the 
mechanical  response  of  the  multi-ply  fabric-based  flexible  armor  material  during  its  impact  with  various 
projectiles  without  significantly  sacrificing  the  key  physical  aspects  of  the  fabric  microstructure,  architec¬ 
ture  and  behavior.  To  validate  the  new  model,  a  comparative  finite  element  method  ( FEM )  analysis  is  car¬ 
ried  out  in  which:  (a)  the  plain-woven  single-ply  fabric  is  modeled  using  conventional  shell  elements  and 
weaving  is  done  in  an  explicit  manner  by  snaking  the  yarns  through  the  fabric  and  (b)  the  fabric  is  treated 
as  a  planar  continuum  surface  composed  of  conventional  shell  elements  to  which  the  new  meso-scale 
unit-cell  based  material  model  is  assigned.  The  results  obtained  show  that  the  material  model  provides 
a  reasonably  good  description  for  the  fabric  deformation  and  fracture  behavior  under  different  combina¬ 
tions  of  fixed  and  free  boundary  conditions. 
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1.  Introduction 

To  respond  to  the  new  enemy  threats  and  warfare  tactics,  mil¬ 
itary  systems,  in  particular  those  supporting  the  US  ground  forces, 
are  being  continuously  transformed  to  become  faster,  more  agile, 
and  more  mobile  so  that  they  can  be  quickly  transported  to  oper¬ 
ations  located  throughout  the  world.  Consequently,  an  increased 
emphasis  is  being  placed  on  the  development  of  improved  light¬ 
weight  body-armor  and  lightweight  vehicle-armor  systems  as  well 
on  the  development  of  new  high-performance  armor  materials. 
High-performance  fiber-based  materials  have  been  exploited  for 
both  body-armor  (e.g.  as  soft,  flexible  fiber  mats  for  personal-ar¬ 
mor  vests)  and  for  the  vehicle-armor  systems  (e.g.  as  reinforce¬ 
ments  in  rigid  polymer  matrix  composites,  PMCs,  for  lightweight 
vehicle-armor  systems). 

Throughout  history,  flexible  lightweight  materials  have  been 
used  in  body-armor  systems  to  provide  protection  against  speci¬ 
fied  threats,  at  reduced  weight  and  without  compromising  person’s 
mobility.  Early  materials  used  included  leather,  silk,  metal  chain 
mail  and  metal  plates.  Replacement  of  metal  with  a  nylon  ( poly¬ 
amide )  fabric  and  an  E-glass  fiber/ethyl  cellulose  composite  in 
body-armor  systems  can  be  linked  to  the  Korean  War  [1]. 
Although,  primarily  due  to  their  low  cost,  nylon  and  E-glass  fibers 
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are  still  being  used  today,  high-performance  polymeric  fibers  (typ¬ 
ically  used  in  the  form  of  woven  fabrics)  are  now  the  standard  in 
most  fiber-reinforced  body-armor  applications.  To  increase  the 
ballistic  performance  of  the  body-armor  vests  relative  to  up  to 
0.30  caliber  threats,  ceramic  insert  strike-plates  are  commonly 
used  [2]. 

High-performance  polymeric  fibers  used  today  are  character¬ 
ized  by  substantially  improved  strength,  stiffness  and  ballistic  per¬ 
formance.  Among  these  high-performance  fibers  the  most  notable 
are:  (a)  poly-aramids  (e.g.  Kevlar®,  Twaron®,  Technora®);  (b)  highly 
oriented  poly-ethylene  (e.g.  Spectra®,  Dyneema®);  (c)  poly-benz- 
obis-oxazole,  PBO  (e.g.  Zylon®),  and  (d)  poly-pyridobisimi-dazolet 
PIPD  (e.g.  M5®).  When  tested  in  tension,  all  these  materials  differ 
significantly  from  the  nylon  fibers,  having  very  high  absolute  stiff¬ 
ness,  extremely  high  density-normalized  strength,  and  quite  low 
(<4%)  strains-to-failure.  These  fibers  essentially  behave,  in  tension, 
as  rate-independent  linear-elastic  materials.  When  tested  in  trans¬ 
verse  compression,  however,  these  fibers  are  similar  to  nylon  and 
can  undergo  large  plastic  deformation  without  a  significant  loss 
in  their  tensile  load-carrying  capacity.  This  behavior  is  quite  differ¬ 
ent  from  that  found  in  carbon  or  glass  fibers,  which  tend  to  shatter 
under  transverse  compression  loading  conditions.  Among  the 
polymeric  fibers  mentioned  above,  the  best  performer,  PBO,  has  a 
tensile-strength  of  ca.  5  GPa  (i.e.  more  than  three  times  the 
tensile-strength  level  of  the  strongest  “ piano-wire ”  steel  but  at  a 
one-fifth  of  the  steel  density,  resulting  in  an  ~1 5-fold  higher 
density-normalized  strength). 
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As  mentioned  above,  fabrics  based  on  high-performance  fibers 
are  extensively  employed  in  variety  of  ballistic  and  impact  protec¬ 
tion  applications.  Despite  the  fact  that  over  the  past  two  decades, 
there  has  been  a  great  deal  of  work  done  on  understanding  the  bal¬ 
listic  behavior  of  these  fabrics  using  various  analytical  and  numer¬ 
ical  techniques,  the  design  of  fabric  armor  systems  remains  largely 
based  on  the  employment  of  extensive  experimental  test  pro¬ 
grams,  empiricism  and  old  practices.  While  such  experimental  pro¬ 
grams  are  critical  for  ensuring  the  utility  and  effectiveness  of  the 
armor  systems,  they  are  generally  expensive,  time-consuming 
and  involve  destructive  testing.  Consequently,  there  is  a  continuing 
effort  to  reduce  the  extent  of  these  experimental  test  programs  by 
complementing  them  with  the  corresponding  computation-based 
engineering  analyses  and  simulations. 

Among  the  main  computational  engineering  analyses  used  to 
model  ballistic  performance  of  flexible  armor  the  following  main 
classes  can  be  identified: 

(a)  Finite  element  analyses  based  on  the  use  of  pin-jointed 
orthogonal  bars  to  represent  flexible-fabric  yarns.  The  most 
notable  studies  falling  into  this  category  of  analyses  are 
those  performed  by  Roylance  and  Wang  [3],  Shim  et  al.  [4], 
Lim  et  al.  [5],  Shahkarami  et  al.  [6],  Johnson  et  al.  [7],  and  Bil¬ 
lon  and  Robinson  [8].  While  the  pin-jointed  orthogonal  bars 
based  finite  element  analyses  have  proven  to  be  very  effi¬ 
cient  in  approximating  the  dynamic  behavior  of  woven  fab¬ 
rics,  the  discrete  nature  of  the  yarn  models  was  associated 
with  inherent  oversimplifications  that  significantly  limited 
the  predictive  capability  of  the  analyses.  In  particular, 
important  contributions  associated  with  the  weave  architec¬ 
ture,  surface-finish  and  friction  governed  yarn-to-yarn  and 
layer-to-layer  contacts  (in  multi-layer  fabrics)  could  not  be 
accounted  for. 

(b)  More-detailed  full-blown  3D  continuum  finite  element  anal¬ 
yses  such  as  the  ones  carried  out  by  Shockey  et  al.  [9],  Duan 
et  al.  [10-13],  Zhang  et  al.  [14],  etc.  have  also  been  investi¬ 
gated.  While  these  analyses  have  proven  to  be  powerful 
tools  for  capturing  and  elucidating  the  detailed  dynamic 
response  of  single-layer  fabrics,  they  are  computationally 
very  demanding  when  applied  to  practical  armor  systems 
which  typically  contain  30-50  fabric  layers/plies. 

(c)  Unit-cell  based  approaches  have  been  used  extensively  in 
order  to  derive  the  equivalent  (smeared)  continuum-level 
(membrane/shell)  material  models  of  textile  composites 
from  the  knowledge  of  the  meso-scale  fiber  and  yarn  prop¬ 
erties,  fabric  architecture  and  inter-yarn  and  inter-ply  fric¬ 
tional  characteristics.  Among  the  most  notable  studies 
based  on  these  analyses  are  those  carried  out  by  Kawabata 
et  al.  [15-17]  who  introduced  simple  analytical  models  to 
capture  the  uniaxial,  biaxial  and  shear  behavior  of  fabrics. 
Furthermore,  Ivanov  and  Tabiei  [18]  proposed  a  micro¬ 
mechanical  material  model  for  a  woven  fabric  (in  which  a 
visco-elastic  constitutive  model  was  used  to  represent  the 
mechanical  behavior  of  the  yarns)  for  the  use  in  non-linear 
finite  element  impact  simulations.  In  deriving  the  material 
model,  Ivanov  and  Tabiei  [18]  considered  the  motion  of 
the  yarn  crossover  point  and  developed  a  procedure  for 
determining  the  equilibrium  position  of  this  point  under 
the  applied  unit-cell  strains.  Recently,  King  et  al.  [19]  pro¬ 
posed  a  new  approach  for  deriving  the  continuum-level 
material  model  for  fabrics  based  on  the  properties  of  the 
yarns  and  the  weave  architecture  which  involves  the  use 
of  an  energy  minimization  technique  to  establish  the  rela¬ 
tionship  between  the  configurations  of  the  fabric  structure 
to  the  microscopic  deformation  of  fabric  components.  Simi¬ 
lar  unit-cell  based  continuum-level  membrane/shell  mate¬ 


rial  models  have  been  developed  by  Boisse  et  al.  [20]  and 
Peng  and  Cao  [21].  Also,  Shahkarami  and  Vaziri  [22]  pro¬ 
posed  a  similar  but  simpler  model  to  that  introduced  by  King 
et  al.  [19]  and  provided  a  detailed  account  of  its  incorpora¬ 
tion  into  a  material  model  subroutine  which  can  be  readily 
coupled  with  commercial  dynamic-explicit  finite  element 
codes. 

(d)  The  use  of  higher-order  membrane/shell  finite  element  anal¬ 
yses  to  represent  the  dynamic  response  of  fabric  under  bal¬ 
listic  loading  conditions  and  overcome  the  aforementioned 
computational  cost  associated  with  the  use  of  full  3D  finite 
element  analyses  of  the  yarn/fabric  structure.  Among  the 
studies  falling  into  this  category,  the  most  notable  is  the 
one  carried  out  by  Scott  and  Yen  [23].  While  the  use  of 
higher-order  membrane  elements  was  found  to  be  indeed 
advantageous  computationally,  it  was  never  fully  validated 
by  comparing  its  results  against  either  those  obtained  exper¬ 
imentally  or  those  obtained  using  full  3D  finite  element 
analyses. 

As  pointed  out  above,  while  major  efforts  have  been  made  in  re¬ 
cent  years  to  develop  sophisticated  numerical  models  capable  of 
elucidating  the  ballistic  performance  of  fabric  armors,  most  of 
these  models  either  lack  computational  efficiency  or  fail  to  capture 
many  physical  aspects  of  the  yarn  and  fabric  architecture  and/or 
contact  dynamic  phenomena.  Hence,  the  main  objective  of  the 
present  work  is  to  develop  an  efficient  shell-based  meso-scale 
mechanics  unit-cell  based  model  that  captures  the  essential  dy¬ 
namic/ballistic  behavior  of  plain-woven  fabric  under  impact-load¬ 
ing  conditions.  The  term  “ meso-scale ”  is  used  to  denote  yarn-level 
millimeter  length  scale  details  of  the  fabric  microstructure/archi¬ 
tecture.  In  other  words,  finer-scale  molecular-level  and  fiber-level 
material  details  are  not  considered  explicitly  and  instead  only  their 
lumped  contributions  are  taken  into  account.  The  “unit-cell”  term 
is  used  to  denote  the  basic  structural  unit  in  a  woven  single-ply 
fabric  so  that  a  fabric  patch  can  be  considered  as  an  in-plane 
assembly  of  such  units.  The  material  model  developed  in  the  pres¬ 
ent  work  is  essentially  an  extension  of  the  model  recently  pro¬ 
posed  by  Shahkarami  and  Vaziri  [22]  and  includes  all  the  cases  of 
in-plane  unit-cell  deformation  modes  and,  for  each  of  these  modes, 
all  the  contact/no-contact  cases  of  the  crossing  yarns.  In  addition, 
three-dimensional  single  unit-cell  finite  element  analyses  of  differ¬ 
ent  in-plane  and  transverse  loading  scenarios  were  carried  out  to 
assess  the  yarn/yarn-contact  force  vs.  over-closure  relationship, 
as  well  as  the  in-plane  and  out-of-plane  shear  moduli.  The  contact 
force  vs.  over-closure  relationship  and  the  in-plane  and  out-of¬ 
plane  shear  unit-cell  moduli  are  then  used  as  input  in  the  meso- 
scale  unit-cell  based  material  model. 

The  organization  of  the  paper  is  as  follows:  details  regarding  the 
computational  procedures  employed  to  develop  a  new  meso-scale 
unit-cell  based  material  model  for  a  prototypical  plain-woven  sin¬ 
gle-ply  fabric  and  the  implementation  of  this  model  into  a  material 
user  subroutine  suitable  for  use  in  commercial  finite  element  pro¬ 
grams  are  presented  in  Section  2.  The  formulation  of  a  simple  pro¬ 
jectile-armor  impact  problem  used  to  validate  the  new  material 
model  is  described  in  Section  3.  Main  results  obtained  in  the  cur¬ 
rent  work  are  presented  and  discuss  in  Section  4.  The  main  sum¬ 
mary  points  and  conclusions  resulting  from  the  present  work  are 
listed  in  Section  5. 


2.  Computational  procedure 

In  this  section,  a  simple  meso-scale  unit-cell  based  material 
model  for  a  plain-woven  single-ply  fabric  is  developed  and  imple¬ 
mented  into  a  material  user  subroutine  suitable  for  use  in  commer- 
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Fig.  1.  Geometrical  representation  of  a  single  unit-cell  for  a  plain-woven  single-ply 
fabric:  (a)  before  and  (b)  after  application  of  a  normal  biaxial  in-plane  deformation. 

cial  finite  element  packages.  A  simple  schematic  of  the  unit-cell 
which  is  used  to  represent  the  plain-woven  single-ply  fabric  struc¬ 
ture/architecture  allotted  to  a  single-yarn  crossover  in  its  initial 
(un-deformed)  configuration  is  depicted  in  Fig.  la.  It  should  be 
noted  that  the  schematic  displayed  in  Fig.  la  is  a  simplification 
of  the  corresponding  unit-cell  model  based  on  the  three-dimen¬ 
sional  fabric  structure/architecture  whose  solid  8-node  finite  ele¬ 


ment  discretization  is  depicted  in  Fig.  2a.  A  comparison  of  the 
models  depicted  in  Figs,  la  and  2a  shows  that  the  two  yarns  are 
simplified  in  Fig.  la  as  two  two-member  truss  elements  each  with 
a  constant  elliptical  cross-section.  A  deformed  configuration  of  the 
simplified  unit-cell  model  is  depicted  in  Fig.  lb. 

Due  to  the  truss-character  of  the  yarn  elements,  the  contact  be¬ 
tween  the  yarns  of  the  crossover  is  reduced  to  a  point  contact.  It 
should  be  noted  that  due  to  the  finite  thickness  of  the  yarns,  the 
yarn/yarn  point  contact  corresponds  to  an  initial-non-zero, 
Hoi  +  H02,  distance  between  the  summit  points  of  the  two  two- 
member  truss  elements,  where  H0i  and  H02  are  the  initial  out-of¬ 
plane  distances  of  the  two  summit  points.  Furthermore,  it  is  as¬ 
sumed  that,  during  deformation,  the  in-plane  location  of  the  con¬ 
tact  point  (i.e.  the  projection  of  the  summit  points  to  the  shell 
surface)  remains  constant  but  the  summit  points  can  move  in  the 
out-of-plane  direction. 

The  initial  lengths  of  the  corresponding  single  truss  members 
are  L0 1  and  L02.  When  the  unit-cell  is  subjected  to  in-plane  biaxial 
normal  strains,  dimensions  of  the  unit-cell  are  changed  from  their 
initial  values,  y0i  and  y02  to  yi(Hi)  and  y2(H2),  respectively.  When 
the  unit-cell  is  stretched  in  a  particular  in-plane  direction,  the  cor¬ 
responding  yarn  is  either  de-crimped/straightened  or  stretched.  In 
the  former  case,  no  tension  is  built  within  the  yarn,  while  in  the 
latter  case  tension  is  created  within  the  yarn  and,  if  sufficiently 
high,  can  cause  yarn  failure.  Also,  the  extent  of  tension  in  a  given 
yarn  is  affected  by  the  failure  status  of  the  other  crossing  yarn.  If 
the  crossing  yarn  is  not  broken  and  is  in  contact  with  the  stretched 
yarn,  tension  will  develop  in  the  two  yarns,  the  extent  of  which  is 
dependent  of  the  extent  of  contact  over-closure  (i.e.  the  reduction 
in  distance  of  the  two  summit  points  relative  to  the  initial  dis¬ 
tance)  and  the  functional  relationship  which  relates  the  contact 
over-closure  with  the  magnitude  of  the  contact  force.  Due  to  their 
low  buckling  resistance,  compression  cannot  be  built  within  the 
yarns  when  the  cell  is  subjected  to  in-plane  normal  compressive 
strains. 


Fig.  2.  Three-dimensional  finite  element  representation  of  the  unit-cell  in  a  plain-woven  single-ply  fabric:  (a)  initial  configuration;  (b)  after  in-plane  biaxial  stretching; 
(c)  after  in-plane  shear;  (d)  after  transverse  shear. 
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2 A.  A  meso-scale  unit-cell  based  material  model  for  plain-woven 
single-ply  fabric 

As  explained  earlier  the  main  objective  of  the  present  work  is  to 
develop  a  meso-scale  unit-cell  based  material  model  for  a  proto¬ 
typical  plain-woven  single-ply  Kevlar  129  fabric  which  can  be  used 
in  large-scale  computational  analyses  of  multi-ply  flexible  armor 
systems.  The  development  of  this  model  is  presented  in  this  sec¬ 
tion.  Since  the  model  is  intended  for  the  use  in  conjunction  with 
shell  finite  elements,  its  development  is  divided  in  three  parts: 
(a)  the  normal  biaxial  in-plane  response  of  the  unit-cell;  (b)  the 
in-plane  shear  response  of  the  cell;  and  (c)  the  out-of-plane  (trans¬ 
verse)  shear  response  of  the  shell.  In  the  following  three  sections, 
the  development  of  each  if  these  three  components  of  the  material 
model  is  presented.  This  is  followed  by  a  brief  discussion  of  the 
implementation  of  the  material  model  in  a  user  material  subrou¬ 
tine  which  is  coupled  with  the  commercial  finite  element  program 
ABAQUS/Explicit  [24]. 

As  mentioned  above,  when  the  two  yarns  are  not  broken,  in 
contact  and  the  distance  between  the  two  summit  points  is  smaller 
than  the  initial  distance,  contact  force  is  developed  between  the 


two  yarns  and  the  extent  of  this  force  depends  on  the  yarn-contact 
over-closure  and  a  functional  relationship  which  relates  the  con¬ 
tact  force  magnitude  to  the  contact  over-closure.  To  determine 
the  functional  relationship  between  the  yarn-contact  over-closure 
dc\  +  dC2  =  dc  and  the  resulting  contact  force  Fc i  a  three-dimen¬ 
sional  solid  FEM  model  of  a  single  unit-cell  is  constructed,  Fig.  2a 
and  subjected  to  the  balanced  biaxial  in-plane  tensile  loads  and 
the  total  surface  contact  force  monitored  is  a  function  of  the 
yarn-contact  over-closure.  An  example  of  deformed  configuration 
of  the  unit-cell  is  shown  in  Fig.  2b  while  the  corresponding  contact 
force  vs.  yarn  over-closure  results  obtained  in  this  analysis  are  dis¬ 
played  in  Fig.  3a.  Following  Shahkarami  and  Vaziri  [22],  the  Fc  vs.  dc 
relation  is  described  using  the  following  function: 


Eq.  (1)  ensures  that  the  contact  force  initially  increases  linearly 
with  the  over-closure  and  as  the  over-closure  begins  to  approach 
the  yarn  thickness  the  contact  force  becomes  excessively  high. 
Using  a  simple  unconstrained  optimization  procedure,  the  three- 
dimensional  solid  FEM  Fc  vs.  dc  results  displayed  in  Fig.  3a  are  fitted 


Out-of-plane  Shear  Strain 


Fig.  3.  Typical  results  obtained  in  three-dimensional  FEM  analyses  of  the  deformation  of  the  unit-cell:  (a)  under  in-plane  biaxial  tension;  (b)  in-plane  shear;  (c)  out-of-plane/ 
transverse  shear. 
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to  the  function  defined  by  Eq.  (1)  and  the  two  parameters  are 
determined  as:  a  =  0.12389  mm  and  b  =  983.39  mN  in  the  absence 
of  yarn/yarn  friction  and  a  =  0.12884  mm  and  b  =  1298.1  mN  in  the 
case  of  yarn/yarn  frictional  coefficient  of  0.5  and  these  values  for  a 
and  b  are  used  in  the  remainder  of  the  paper.  It  should  be  noted 
that  the  data  pertaining  to  the  case  of  yarn/yarn  frictional  coeffi¬ 
cient  of  0.5  in  Fig.  3a  are  shifted  upward  by  500  mN  to  improve 
clarity  of  the  figure. 

2  A  A.  In-plane  biaxial  behavior  of  the  unit-cell 

When  the  unit-cell  is  subjected  to  in-plane  biaxial  strains,  sev¬ 
eral  cases  have  to  be  considered  depending  on  whether  the  strains 
are  tensile  or  compressive  and  whether  any  of  the  yarns  is  broken: 

(a)  When  both  yarns  are  broken  within  a  cell,  the  corresponding 
finite  element  is  eroded,  i.e.  it  is  removed  from  the  analysis. 

(b)  When  one  of  the  yarns  is  broken,  then  the  other  yarn  is  trea¬ 
ted  as  being  the  only  yarn  in  the  unit-cell.  If  the  un-broken 
yarn  is  subjected  to  compression  or  the  unit-cell  length  in 
the  direction  of  this  yarn,  y{  (z  =  1  or  2  and  pertains  to  the 
un-broken  yarn)  is  smaller  than  the  corresponding  initial 
yarn  length  L0l,  then  no  load  is  supported  by  this  yarn.  In 
other  words,  zero  in-plane  normal  stresses  exist  along  the 
axis  of  both  yarns.  When  yt  >  Loi<  then  tension  is  developed 
within  the  un-broken  yarn.  This  tension  is  defined  by  the  fol¬ 
lowing  single-yarn  expression  as: 

Ti  =  EAi(Li-Loi)/Loi  (2) 

where  E  is  the  yarn  axial  Young’s  modulus,  is  the  cross-sec¬ 
tional  area,  L,  =  y,-  and  i  refers  to  the  un-broken  yarn. 

(c)  When  both  yarns  are  intact,  their  respective  unit-cell  in¬ 
plane  normal  strains  are  both  positive  and  there  is  a  yarn- 
contact  over-closure,  then  the  effect  of  the  contact  force 
has  to  be  taken  into  account  before  the  yarn  tensions  (and 
in  turn,  the  corresponding  in-plane  unit-cell  normal  stres¬ 
ses)  can  be  computed.  A  simple  procedure  for  the  calculation 
of  the  yarn  tensions  and  unit-cell  in-plane  normal  stresses  is 
provided  below.  It  should  be  noted  that  the  total  yarn-con- 
tact  over-closure  is  defined  as,  i.e.  dc  =  dc i  +  dc2.  Using  single 
geometrical  arguments  and  Fig.  la  and  b,  the  current  yarn 
crimp  heights  can  be  defined  as: 

Hi(dci)  =  Hoi  -  dci  (i  =1,2)  (3) 

Likewise,  the  current  single-member  truss  element  lengths 
can  be  defined  as: 

mi)  =  +  (i  =1,2)  (4) 

The  corresponding  yarn  tensions  are  then  defined  by  Eq.  (2). 
Using  the  force  equilibrium  condition  at  the  yarn  summit 
points,  the  following  expression  can  be  derived  for  the  corre¬ 
sponding  contact  forces: 

Fci(dci )  =  2 Tj  ( dci ) Hi  (da)/ Lt  ( dci )  (i  =  1,2)  (5) 

The  yarn/yarn  forces  at  the  two  summit  points,  Fci  and  Fc2, 
defined  by  Eq.  (5)  are  equal  to  each  other  and  are  also  equal 
to  the  contact  force,  Fc,  defined  by  Eq.  (1 ).  Thus  the  following 
system  of  two  non-linear  algebraic  equations  with  two  un¬ 
knowns,  dci  and  dc2,  can  be  defined  as: 

2T\  (dcl )  (Hi  (dcl )  /  L\  (dcl ))  =  2  T2(dc2)(H2(dc2)/L2(dc2))  (6a) 
2T1(dcl)(H1(dcl)/L1(dcl))  =  b(dcl  +  dc2)/(a  -  dcl  -  dc2)  (6b) 

where  r,(dd),  Lj(dd)  and  Hj(dd)  are  defined  by  Eqs.  (2)-(4), 
respectively.  Eqs.  (6a)  and  (6b)  can  be  readily  solved  using 


a  Newton-Raphson  routine.  Once  dc\  and  dc2  are  determined, 
the  two  yarn  tensions  can  be  computed  using  Eq.  (2)  and  in 
turn,  the  two  in-plane  unit-cell  normal  stresses  can  be  deter¬ 
mined  as: 

(ti^Tt/Ai  (i  =  1,2)  (7) 

(d)  When  both  yarns  are  intact  and  one  of  the  in-plane  unit-cell 
normal  strains  is  tensile  while  the  other  one  is  compressive, 
two  possible  cases  can  arise:  (i)  if  the  current  summit-point 
distance  is  larger  than  its  initial  value  (H0i  +  H02 ),  then  the 
yarns  are  not  in  contact  and  the  yarn  which  is  aligned  in 
the  compression-strain  direction  (the  “ compressed  yarn") 
does  not  carry  any  load  and  the  normal  stress  in  the  com¬ 
pression-strain  direction  is  zero.  As  far  as  the  other  direction 
is  concerned,  tension  can  be  developed  in  the  yarn  (when 
yz-  >  L0l,  where  i  pertains  to  the  yarn  aligned  in  the  tensile- 
strain  direction  of  the  unit-cell,  the  “stretched  yarn").  In  this 
case,  the  yarn  tension  is  given  by  the  single-yarn  tension 
relation,  Eq.  (2).  Otherwise,  no  tension  is  built  within  the 
yarn  and  the  corresponding  in-plane  normal  unit-cell  stress 
also  becomes  zero;  (ii)  when  yarn-contact  over-closure  is 
present,  then  tension  builds  up  in  both  yarns  and  the  two 
in-plane  unit-cell  normal  stresses  are  both  positive/tensile. 
To  determine  the  tensions  in  the  two  yarns,  a  similar  proce¬ 
dure  can  be  developed  as  the  one  used  under  point  (c)  in  this 
section.  The  main  difference  is  that  the  reference  configura¬ 
tions  from  which  the  development  of  the  yarn/yarn-contact 
forces  is  measured  are  different  in  the  two  cases.  In  the  case 
of  biaxial  stretching  of  the  unit-cell,  the  reference  configura¬ 
tion  corresponds  to  the  initial  configuration  of  the  fabric  in 
which  the  two  yarns  are  in  contact  but  no  over-closure 
(i.e.  the  contact  force)  is  present.  In  the  present  case,  it  is 
advantageous  to  temporarily  decouple  the  behavior  of  the 
two  yarns.  The  summit-point  height  of  the  “compressed” 
yarn  has  increased  under  unit-cell  compression.  The  sum¬ 
mit-point  of  the  “stretched”  yarn  is  decreased  due  to  unit¬ 
cell  tension  and  the  magnitude  of  this  decrease  is  larger  than 
the  summit-point  height  increase  of  the  “compressed”  yarn. 
Hence,  the  reference  configuration  in  this  case  corresponds 
to  the  one  in  which  the  decrease  in  the  summit-point  height 
of  the  “stretched”  yarn  is  temporarily  set  equal  to  the  sum¬ 
mit-point  height  increase  of  the  “compressed”  yarn.  Then, 
further  decrease  in  the  summit-point  height  of  the 
“stretched”  yarn  will  cause  a  further  increase  in  the  sum¬ 
mit-point  height  of  the  “compressed”  yarn  which  leads  to 
a  contact  force  build-up  and  to  the  development  of  tension 
in  both  yarns.  Once  account  is  taken  of  the  new  reference 
configuration,  yarn  tensions  and  the  corresponding  unit-cell 
normal  stresses  can  be  calculated  using  the  relations  identi¬ 
cal  to  those  presented  under  point  (c)  in  this  section. 


2 A. 2.  In-plane  shear  behavior  of  the  unit-cell 

The  meso-scale  material  model  presented  in  the  previous  sec¬ 
tion  provides  a  way  for  assessing  in-plane  normal  stresses  of  the 
unit-cell.  Since  no  friction  between  the  two  two-member  truss  ele¬ 
ments  is  taken  into  account  in  Fig.  1,  the  unit-cell  model  could  not 
be  used  to  assess  its  in-plane  shear  behavior.  To  overcome  this  lim¬ 
itation,  the  three-dimensional  unit-cell  finite  element  model  dis¬ 
played  in  Fig.  2a  is  subjected  to  an  increase  in-plane  shear  and 
the  corresponding  shear  stress  is  determined  by  dividing  the  total 
in-plane  shear  force  by  the  corresponding  yarn  cross-sectional  area. 
An  example  of  deformed  configuration  of  the  unit-cell  is  displayed 
in  Fig.  2c  while  the  corresponding  shear  stress  vs.  shear  strain  re¬ 
sults  obtained  in  this  analysis  are  displayed  in  Fig.  3b.  The  FEM  re- 
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suits  displayed  in  Fig.  3b  are  next  fitted  to  an  inverse  tangent  func¬ 
tion,  denoted  by  a  solid  line  in  Fig.  3b  to  obtain  a  relationship  be¬ 
tween  the  in-plane  shear  modulus,  Gin.piane ,  and  the  in-plane 
shear  strain,  sin.piane •  The  following  functional  relations  were 
obtained:  Gin_p/ane  (GPa)  =  38.45  +  23.744a  tan(67.0£in_p/ane  -  26.8) 
in  the  absence  of  yarn/yarn  friction,  and  Gin_piane  (GPa)  = 
41.142  +  22.26a  tan(69.0e,n-p/ane  -  27.6)  in  the  case  of  yarn/yarn 
frictional  coefficient  of  0.5.  It  should  be  noted  that  the  data  pertain¬ 
ing  to  the  case  of  yarn/yarn  frictional  coefficient  of  0.5  in  Fig.  3b  are 
shifted  upward  by  10  GPa  to  improve  clarity  of  the  figure. 

2.1.3.  Out-of-plane  (transverse  shear)  behavior  of  the  unit-cell 

Following  a  procedure  similar  to  the  one  described  in  the  previ¬ 
ous  section,  the  three-dimensional  FEM  model  of  the  unit-cell  dis¬ 
played  in  Fig.  2a  is  subjected  to  an  out-of-plane  shear  and  the 
corresponding  shear  stress  determined  by  dividing  the  applied 
transverse  force  by  the  yarn  cross-sectional  area.  An  example  of 
deformed  configuration  of  the  unit-cell  is  displayed  in  Fig.  2d  while 
the  corresponding  transverse  shear  modulus  vs.  transverse  shear 
strain  results  obtained  in  this  analysis  are  displayed  in  Fig.  3c. 
The  results  displayed  in  Fig.  3c  are  fitted  to  a  constant  relation 
yielding  the  transverse  shear  moduli  of  0.89  and  0.97  GPa  for  the 
cases  of  yarn/yarn  frictional  coefficient  of  0.0  and  0.5,  respectively. 
It  should  be  noted  that  the  data  pertaining  to  the  case  of  yarn/yarn 
frictional  coefficient  of  0.5  in  Fig.  3c  are  shifted  upward  by  0.1  GPa 
to  improve  clarity  of  the  figure. 

2.2.  Material  model  implementation  in  a  User-material  Subroutine 

The  meso-scale  unit-cell  based  material  model  described  in  the 
previous  section  is  next  implemented  in  the  material  user  subrou¬ 
tine,  VUMAT,  of  the  commercial  finite  element  program  ABAQUS / 
Explicit  [24].  This  subroutine  is  compiled  and  linked  with  the  finite 
element  solver  and  enables  ABAQUS/Explicit  to  obtain  the  needed 
information  regarding  the  state  of  the  material  and  the  material 
mechanical  response  during  each  time  step,  for  each  integration 
point  of  each  element.  In  the  present  work  first-order  4-node  gen¬ 
eral-purpose  reduced-integration  shell  elements  (ABAQUS/Explicit 
designation  S4R)  are  used.  Due  to  the  use  of  the  Simpson’s  numer¬ 
ical-integration  method  for  the  calculation  of  through-the-thick- 
ness  deformation  response  of  the  shell,  an  odd  number  (greater 
than  1)  of  integration  points  has  to  be  used  in  the  through-the- 
thickness  direction.  The  results  obtained  in  the  present  work  sug¬ 
gest  that  selecting  three  integration  points  provides  a  good  compro¬ 
mise  between  computational  efficiency  and  accuracy. 

The  essential  features  of  the  coupling  between  the  ABAQUS/Ex¬ 
plicit  finite  element  solver  and  the  VUMAT  Material  User  Subrou¬ 
tine  at  each  time  increment  at  each  integration  point  of  each 
element  can  be  summarized  as  follows: 

(a)  The  corresponding  previous  time  increment  stresses  and 
material  state  variables  as  well  as  the  current  time  step 
strain  increments  are  provided  by  the  ABAQUS/Explicit  finite 
element  solver  to  the  material  subroutine.  In  the  present 
work,  the  unit-cell  strains,  warp  and  weft  yarn  summit- 
point  heights,  yarn  lengths,  as  well  as  the  flags  pertaining 
to  the  failure  status  of  the  yarns  and  the  deletion  status  of 
the  element  are  used  as  state  variables. 

(b)  Using  the  information  provided  in  (a),  and  the  meso-scale 
unit-cell  used  material  model  presented  in  the  previous  sec¬ 
tion,  the  material  stress  state  as  well  as  values  of  the  mate¬ 
rial  state  variables  at  the  end  of  the  time  increment  are 
determined  within  the  VUMAT  and  returned  to  the  ABA¬ 
QUS/Explicit  finite  element  solver.  In  addition,  the  changes 
in  the  total  internal  and  the  inelastic  energies  (where  appro¬ 
priate)  are  computed  and  returned  to  the  solver. 


It  should  be  noted  that  within  the  VUMAT  only  the  normal  and 
shear  in-plane  response  of  the  shell  material  is  computed.  Trans¬ 
verse  shear  stiffness  of  the  shell  elements,  as  computed  using  the 
three-dimensional  FEM  procedure  outlined  in  the  previous  section, 
has  to  be  defined  as  part  of  the  overall  FEM  model  definition  out¬ 
side  the  VUMAT.  Using  the  provided  values  for  the  transverse  shear 
modulus,  ABAQUS/Explicit  uses  a  simple  procedure  in  which  the 
transverse  response  of  a  shell  is  approximated  by  the  transverse 
response  of  an  analogous  solid  finite  element. 

To  provide  a  more  detailed  insight  into  the  implementation  of 
the  current  meso-scale  unit-cell  based  material  model  in  the  VU¬ 
MAT,  a  flow  chart  is  provided  in  Fig.  4.  It  is  important  to  note  that 
ABAQUS/Explicit  provides  tensorial  shear  strain  rather  than  the 
engineering  shear  strain  increments.  The  procedure  used  to  com¬ 
pute  unit-cell  stresses  at  the  end  of  the  current  time  step  is  clearly 
dependent  on  failure  status  of  the  two  yarns  at  the  end  of  the  pre¬ 
vious  time  step  and  on  the  current  (tension  vs.  compression)  state 
of  the  unit-cell  edges.  Yarns  cannot  support  compressive  loads  and 
broken  yarns  cannot  support  tensile  loads.  In  both  cases  the  asso¬ 
ciated  in-plane  normal  stress  components  are  zero. 

3.  Validation  of  the  material  model 

As  stated  earlier,  the  main  objective  of  the  present  work  was  to 
develop  a  fabric-structure/architecture  dependent  and  physically 
based,  computationally  efficient  material  model  which  can  be  used 
in  ballistic-protection  performance  analyses  of  multi-ply  armor.  In 
this  section,  a  simple  projectile/armor  impact  problem  is  de¬ 
scribed.  The  problem  is  used  to  carry  out  preliminary  testing  and 
validation  of  the  meso-scale  unit-cell  based  material  model  for 
plain-woven  fabric  developed  in  the  present  work. 

The  initial  configurations  of  the  projectile/armor  finite  element 
systems  analyzed  here  are  shown  as  Fig.  5a  and  b.  In  both  cases,  a 
rigid  spherical  projectile  with  a  4  mm  radius  and  a  2.12  g  weight  is 
propelled  at  an  initial  velocity  of  400  m/s  in  the  direction  normal 
to  the  single-ply  armor  surface  and  toward  the  center-point  of 
the  armor.  The  armor  is  modeled  as  a  36  mm  by  36  mm  square 
patch.  In  Fig.  5a,  the  plain-woven  fabric  structure  is  modeled 
explicitly  by  snaking  through  orthogonally  oriented  warp  and  weft 
yarns.  The  square-shaped  fabric  patch  contained  44  warp  and  44 
weft  yarns.  Yarns  are  considered  to  have  a  constant  0.040  mm2 
hexagonal  cross-sectional  area.  The  hexagonal  cross-section  was 
used  since  this  is  the  best  approximation  to  the  actual  elliptical 
cross-section  which  can  be  obtained  using  two-element  wide 
yarns  and  the  “Nodal  Thickness ”  option  available  in  ABAQUS/Expli¬ 
cit.  The  yarn  width  is  set  to  0.615  mm,  the  peak  height  to 
0.105  mm,  the  crimp  height  to  0.088  mm  and  crimp  wave-length 
to  1.64  mm.  All  these  values  are  consistent  with  those  used  by 
Duan  et  al.  [10-13]  and  Zhang  et  al.  [14]  and,  for  an  effective  yarn 
density  of  0.6  g/cm3,  they  yield  an  effective  fabric  areal-density  of 
3.31  g/cm2. 

The  yarn-material  is  assumed  to  be  linear-elastic  orthotropic  (or 
more  precisely,  transversely  isotropic)  with  the  unique  material 
direction  being  aligned  in  the  yarn-axis  direction).  The  orthotropic 
linear-elastic  yarn-material  properties  used  are  listed  in  Table  1 
and  they  relate  to  the  corresponding  properties  of  Kevlar®  129. 
Low  values  for  the  transverse  normal  and  shear  moduli  and  for 
the  Poisson’s  ratios  listed  in  Table  1  arise  from  the  fact  that  the  fi¬ 
bers  bundled  within  yarns  are  only  weakly  coupled  to  each  other. 

In  the  case  of  Fig.  5b,  the  armor  is  modeled  as  a  constant-thick¬ 
ness  shell-based  continuous  surface.  In  other  words,  warp  and  weft 
yarns  and  their  weaving  are  not  modeled  explicitly.  The  effect  of 
yarn  weaving,  however,  is  included  implicitly  through  the  use  of 
the  meso-scale  unit-cell  based  material  model  which  is  assigned 
to  each  shell  element  of  the  fabric  in  Fig.  5b. 
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Fig.  4.  The  flow  chart  used  during  implementation  of  the  meso-scale  unit-cell  based  material  model  in  a  User-material  Subroutine. 


Fig.  5.  The  initial  configurations  of  the  projectile/fabric  systems  for  (a)  a  shell  FEM  analysis  in  which  yarn  weaving  is  modeled  explicitly;  and  (b)  a  shell  FEM  analysis  in  which 
the  effects  of  yarn  weaving  are  included  implicitly  through  the  use  of  an  meso-scale  unit-cell  based  material  model. 


Three  different  types  of  far-held  boundary  conditions  applied  to 
the  edges  of  the  fabric  patch  were  considered:  (a)  all  four  fabric 
edges  are  clamped;  (b)  two  opposite  fabric  edges  are  clamped 
and  the  other  two  are  free;  (c)  all  four  fabric  edges  are  set  free. 


A  simple  penalty-based  algorithm  is  used  to  model  yarn/yarn 
and  projectile/fabric  interactions.  As  simple  Coulomb  model  was 
used  to  analyze  yarn/yarn  and  projectile/fabric  friction.  Two  fric¬ 
tional  conditions  were  considered:  (a)  both  the  yarn/yarn  friction 
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Table  1 

The  orthotropic  linear-elastic  material  data  (GPa)  for  yam  [11]. 


E  ii 

E22 

E33 

G12 

G13 

G23 

V12 

V13 

V23 

164.0 

3.28 

3.28 

3.28 

3.28 

3.28 

0.0 

0.0 

0.0 

coefficient  Hyjy  and  the  projectile/fabric  friction  coefficient  /ip//  are 
set  to  0.5;  and  (b)  no  yarn/yarn  or  projectile/fabric  friction. 

Due  to  explicit  account  of  the  warp  and  weft  yarns  and  their 
weaving  which  entails  the  use  of  finer  meshes,  the  FEM  model  dis¬ 
played  in  Fig.  5a  is  typically  computationally  3-4  times  more 
expensive  than  the  one  displayed  in  Fig.  5b. 

To  validate  the  meso-scale  unit-cell  based  material  model,  the 
interaction  between  the  projectile  and  the  armor  is  analyzed  using 
the  two  FEM  model  and  two  aspects  of  these  interactions  are  clo¬ 
sely  examined:  (a)  the  residual  velocity  of  the  projectile  after  the 
projectile  has  successfully  penetrated  the  armor  and  (b)  the  extent, 
the  temporal  evolution  and  the  spatial  distribution  of  fabric  armor 
damage. 

It  should  be  noted  that  the  impact  of  a  projectile  with  the  fabric 
is  associated  with  the  initiation  of  several  phenomena,  the  most 
important  of  which  are:  (a)  a  resisting  force  is  exerted  by  the  fabric 
onto  the  projectile  which  causes  a  reduction  in  the  projectile  veloc¬ 
ity;  (b)  at  the  same  time,  the  fabric  is  being  deformed  and  acceler¬ 
ated;  (c)  strain  waves  generated  in  the  impact  region  propagate 
along  the  yarns  toward  the  fabric  edges,  where  they  are  reflected. 
In  the  absence  of  any  external  force  acting  on  the  projectile/fabric 
system,  the  total  energy  of  the  system  must  be  conserved.  In  gen¬ 
eral,  the  energy  dissipation  due  to  projectile  deformation,  fiber 
intermolecular  friction,  wind  resistance  and  acoustic  losses  are 
all  assumed  to  be  negligible.  Consequently,  any  loss  in  projectile 
kinetic  energy  A Epk  can  be  mainly  attributed  to  the  following  three 
energy-absorbing  mechanisms:  yarn  strain  energy  Eys,  yarn  kinetic 
energy  Eyk,  and  the  energy  lost  due  to  frictional-sliding  Ef. 

The  energy  conservation  principle  then  requires  that: 

AEpk  =  Eys  +  Eyk  +  Ef  (8) 

As  clearly  shown  by  Duan  et  al.  [10-13],  the  loss  of  projectile  ki¬ 
netic  energy  A Epk  is  governed  (through  the  three  aforementioned 
energy-absorption  mechanisms)  by  several  factors  such  as  the 
material  properties  of  the  constituent  fibers,  fabric  structure, 
boundary  conditions,  projectile  geometry,  impact  velocity,  friction 
between  the  projectile  and  the  fabric,  and  yarn-to-yarn  and  fiber- 
to-fiber  friction  within  the  fabric  itself. 

Traditionally,  the  overall  decrease  in  projectile  kinetic  energy, 
A  Epk,  is  determined  using  ballistic  impact  experiments  in  which 
only  the  projectile  initial  velocity  and  the  residual  velocity  vr 
are  measured.  The  overall  decrease  in  projectile  kinetic  energy, 
A Epk,  is  then  defined  as: 

A  Epk=^m[vt-v2r]  (9) 

where  m  is  the  mass  of  the  projectile.  Recently,  Starratt  et  al.  [25] 
developed  a  ballistic  impact  method  for  continuous  measurement 
of  projectile  velocity  v(t).  The  loss  of  projectile  kinetic  energy  as  a 
function  of  time  t  during  impact  is  then  defined  as: 

A  Epk=^m[vj-v(t)2}  (10) 

While  this  experimental  technique  may  greatly  help  improve 
our  understanding  of  the  ballistic  impact  behavior  of  fabrics,  many 
essential  physical  phenomena  accompanying  projectile/fabric 
interactions  can  only  be  obtained  using  numerical  analyses  and 
simulations  like  the  ones  used  in  the  present  work. 


4.  Results  and  discussion 

In  this  section,  selected  results  are  presented  and  discussed.  As 
mentioned  earlier,  the  primary  purpose  of  the  present  work  was 
to  develop,  implement  and  validate  a  computationally  efficient, 
meso-scale  unit-cell  based  material  model  for  plain-woven  sin¬ 
gle-ply  fabric  armor.  Since  the  key  functional  requirement  for 
an  armor  system  is  to  absorb  the  kinetic  energy  carried  by  the 
projectile,  a  quantitative  comparison  of  the  results  pertaining  to 
the  temporal  evolution  of  the  absorbed  energies  (through  yarn 
deformation  and  fracture,  fabric  acceleration  and  frictional-sliding 
based  energy  dissipation)  obtained  using  the  two  FEM  formula¬ 
tions  discussed  in  the  previous  section  is  presented  in  the  subse¬ 
quent  sections.  The  results  obtained  for  the  case  of  far-held 
boundary  conditions  corresponding  to  the  all  four  fabric  edges 
being  clamped,  and  for  the  two  frictional  conditions  are  presented 
and  discussed  in  the  next  section.  The  effects  of  the  far-held 
boundary  conditions  are  then  discussed  in  the  subsequent 
section. 

4A.  Four-edges  clamped  far- field  boundary  conditions 

All  the  results  presented  in  this  section  correspond  to  the  case 
of  hxed  boundary  conditions  applied  to  all  four  edges  of  the  fabric. 
As  discussed  earlier,  yarn/yarn  friction  and  projectile/fabric  friction 
are  generally  present  in  experimental  and  held  tests  of  flexible- 
fabric  armor,  and  as  shown  by  Tan  et  al.  [24]  in  a  series  of  post-im¬ 
pact  fabric-inspection  studies,  play  an  important  role  in  absorbing 
the  projectile  kinetic  energy.  For  example,  in  the  presence  of  fric¬ 
tion  yarn  crossover  domains  near  the  impact  region  are  found  to 
be  characterized  by  extensive  hber  breakage,  while  reduced  fric¬ 
tion  promotes  slippage  between  the  warp  and  weft  yarns  and  typ¬ 
ically  the  impact-induced  fabric-perforation  is  smaller  than  the 
projectile  size.  Consequently,  the  hrst  case  analyzed  here  is  the 
one  in  which  both  the  yarn/yarn  friction  coefficient  /iy/y  and  the 
projectile/fabric  friction  coefficient  /rp//  are  set  to  a  non-zero 
(=0.5)  value. 

4AA.  Fabric  deformation  and  yarn  fracture  in  the  presence  of  friction 

It  is  well-established  that  when  a  projectile  hits  the  armor,  two 
waves  are  generated:  (a)  a  longitudinal  wave  which  travels  away 
from  the  point  of  impact  along  the  principal  yarns  (the  yarns  di¬ 
rectly  hit  by  the  projectile)  and  along  those  secondary  yarns  (the 
yarns  not  directly  impacted  by  the  projectile)  which  are  interacting 
with  the  principal  yarns.  Propagation  of  the  longitudinal  wave  out¬ 
ward  from  the  point  of  impact  enables  a  large  fraction  of  the  fabric 
armor  to  undergo  deformation  and,  thus,  absorb  the  kinetic  energy 
of  the  projectile;  and  (b)  a  transverse  wave  which  propagates  out¬ 
ward  from  the  point  of  impact  at  a  velocity  substantially  lower 
than  the  sound  speed  and  is  responsible  for  stretching  of  the  prin¬ 
cipal  and  the  “engaged”  secondary  yarns.  Since  the  sound  speed  is 
controlled  by  the  yarn  axial  stiffness  and  density,  and  these  two 
quantities  are  identical  in  the  two  FEM  models,  the  temporal  and 
spatial  evolutions  of  the  longitudinal  wave  front  are  found  to  be 
quite  comparable  in  the  two  models  (the  results  not  shown  for 
brevity).  The  temporal  evolution  and  the  spatial  distribution  of 
the  transverse  wave  front  are  affected  to  a  greater  extent  by 
yarn/yarn  interactions.  Therefore,  a  comparison  of  the  transverse- 
deflection  wave  front  propagation  results  obtained  using  the  two 
FEM  analyses  can  be  considered  as  a  good  validity  test  for  the 
meso-scale  unit-cell  based  material  model  developed  in  the  pres¬ 
ent  work. 

Examples  of  the  temporal  evolution  of  deformation  within  the 
fabric  obtained  using  the  FEM  analysis  involving  discrete  orthogo¬ 
nal  yarns  plain-woven  into  a  single-ply  fabric  (referred  to  as  the 


3698 


M.  Grujicic  et  al. / Materials  and  Design  30  (2009)  3690-3704 


Fig.  6.  The  temporal  evolution  of  deformation  in  the  fabric  for  the  yarn-level  FEM  model  under  the  yarn/yarn  jiyiy  and  projectile/fabric  ^//friction  coefficients  of  0.5.  Contour 
bands  correspond  to  different  values  of  the  transverse  displacement,  i.e.  the  displacements  normal  to  the  fabric  surface.  All  four  fabric  edges  are  fixed. 


“yarn-level  FEM  analysis ”)  and  the  FEM  analyses  based  on  the 
meso-scale  unit-cell  material  model  (referred  to  as  the  “ unit-cell 
based  FEM  analysis ”)  are  displayed  in  Figs.  6a-d  and  7a-d,  respec¬ 
tively.  In  these  figures,  side  and  top  views  of  the  fabric  along  with 
superimposed  contour  plots  of  the  transverse  displacement  (the 
displacement  normal  to  the  fabric  surface)  are  shown.  A  simple 
comparison  of  the  results  displayed  in  Figs.  6a-d  and  7a-d  reveals 
that  the  temporal  evolution  of  the  deformation  state  of  fabric  is 
quite  similar  in  the  two  analyses  and  can  be  summarized  as 
follows: 

(a)  Initially,  the  shape  of  the  transverse-deflection  wave  front  in 
the  fabric  is  nearly  circular  in  both  analyses  and,  thus,  essen¬ 
tially  identical  to  the  shape  of  the  projectile/fabric  contact- 
zone,  Figs.  6a  and  7a. 

(b)  In  the  case  of  the  yarn-level  FEM  analysis,  as  the  time  pro¬ 
ceeds,  the  transverse-deflection  wave  generated  within  the 


principal  yarns  (the  yarns  which  are  in  direct  contact  with 
the  projectile)  propagates  outward  and,  through  their  inter¬ 
actions  with  the  secondary  yarns  (the  yarns  which  are  not  in 
direct  contact  with  the  projectile),  at  the  yarn  crossovers, 
cause  the  secondary  yarns  to  also  deflect  in  the  transverse 
direction.  Consequently,  the  transverse  wave  front  begins 
to  acquire  a  near  square  shape,  with  the  square  center  coin¬ 
ciding  with  the  impact-zone  center,  Fig.  6b  and  c.  In  the  case 
of  the  unit-cell  based  FEM  analysis,  the  fabric  is  modeled  as  a 
continuous  surface  and,  hence,  the  transverse  wave  can 
readily  propagate  in  each  radial  direction  from  its  source 
point  (i.e.  from  the  point  of  initial  impact).  Consequently, 
the  transverse-deflection  wave  front  acquires  a  nearly  circu¬ 
lar  (more  precisely,  a  polygonal)  shape,  Fig.  7b  and  c. 

(c)  Despite  the  differences  in  the  transverse-deflection  wave 
front  shape  discussed  in  (b),  the  overall  size  of  the  fabric 
region  which  has  undergone  transverse  deflections  and  the 
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Fig.  7.  The  temporal  evolution  of  deformation  in  the  fabric  for  the  unit-cell  based  FEM  model  under  the  yarn/yarn  /ry/y  and  projectile/fabric  /rp//  friction  coefficients  of  0.5. 
Contour  bands  correspond  to  different  values  of  the  transverse  displacement,  i.e.  the  displacements  normal  to  the  fabric  surface.  All  four  fabric  edges  are  fixed. 


extents  of  these  deflections  are  quite  comparable  in  the  two 
FEM  analyses. 

(d)  The  square-shaped  transverse-deflection  wave  front  contin¬ 
ues  to  propagate  towards  the  clamped  edges  of  the  fabric 
and,  upon  reaching  the  edges,  it  is  reflected  back  towards 
the  center  of  the  impact-zone  and  the  wave  front  acquires 
an  octagonal  shape,  Fig.  6d.  The  more  circle-shaped  trans¬ 
verse-deflection  wave  front  undergoes  a  similar  shape 
change  upon  reaching  and  reflecting  from  the  fabric  edges, 
Fig.  7d. 

(e)  The  first  evidence  of  fabric  failure  is  seen  in  the  case  of  the 
unit-cell  based  FEM  analysis,  Fig.  7b.  At  the  same  post¬ 
impact  time,  the  yarns  underneath  the  projectile  in  the  case 
of  the  yarn-level  FEM  analysis  undergo  considerable  lateral 
motion  which  reduces  the  extent  of  tension  within  the  yarns 
and  delays  their  failure,  Fig.  6b.  However,  at  longer  post¬ 
impact  times,  Figs.  6c  and  d  and  7c  and  d,  the  extent  of  fabric 
damage  becomes  quite  comparable  in  the  two  analyses. 


(f)  The  final  fabric  penetration  hole  size  (measured  in  the  yarn 
directions)  is  somewhat  smaller  in  the  case  of  yarn-level 
analysis,  Fig.  6d,  than  in  the  unit-cell  based  case,  Fig.  7d.  This 
relatively  small  discrepancy  is  related  to  the  ability  of  some 
of  the  yarns  which  were  pushed  laterally  by  the  projectile  to 
recoil.  Such  elastic  recovery  was  also  observed  in  the  case  of 
the  continuous  fabric,  but  to  lower  extent.  Nevertheless,  the 
size  of  the  penetration  holes  are  quite  comparable  in  the  two 
cases  and  considering  the  fact  that  the  FEM  mesh  is  signifi¬ 
cantly  finer  and  can  capture  more  details  in  the  case  of  the 
yarn-level  FEM  model,  the  overall  agreement  between  the 
two  analyses  relative  to  the  extent  of  damage  in  the  armor 
can  be  deemed  as  reasonably  good. 

(g)  At  ~25-27  ps,  the  projectile  completely  penetrates  the  fab¬ 
ric  and  continues  to  move  at  a  residual  velocity  of 
393.3  m/s  in  the  case  of  yarn-level  FEM  analysis  and  at  a 
velocity  of  393.7  m/s  in  the  case  of  unit-cell  based  FEM 
analysis. 
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Fig.  8.  The  temporal  evolution  of  deformation  in  the  fabric  for  the  yarn-level  FEM  model  under  the  yarn/yarn  /ry/y  and  projectile/fabric  /rp//  friction  coefficients  of  0.0.  Contour 
bands  correspond  to  different  values  of  the  transverse  displacement,  i.e.  the  displacements  normal  to  the  fabric  surface.  All  four  fabric  edges  are  fixed. 


(h)  It  should  be  noted  that  despite  the  fact  that  the  projectile/ 
fabric  model  has  two  vertical  planes  of  symmetry,  the  dam¬ 
age  region  is  asymmetric  in  both  types  of  analyses.  The  rea¬ 
son  for  this  is  that  the  finite  element  discretization  of  the 
projectile  into  tetrahedron  elements  did  not  possess  two 
planes  of  symmetry.  The  results  displayed  in  Figs.  6a-d 
and  7a-d  thus  also  reveal  the  effect  of  small  geometrical  per¬ 
turbations  in  the  (spherical)  projectile  on  the  fabric  failure 
response. 


4.1.2.  Fabric  deformation  and  yarn  fracture  in  the  absence  of  friction 
The  results  displayed  in  Figs.  6a-d  and  7a-d  and  their  discus¬ 
sion  presented  in  the  previous  section  clearly  establish  that,  at 
the  yarn/yarn  and  projectile/yarn  friction-coefficient  conditions 
of  0.5,  the  two  FEM  analyses  yield  comparable  results  regarding 
the  temporal  evolution  of  the  deformation  within  the  fabric  during 
impact.  However,  some  minor  discrepancies  were  observed  rela¬ 


tive  to  the  extent  of  fabric  damage.  In  this  section,  the  correspond¬ 
ing  computational  results  obtained  using  the  two  FEM  analyses  for 
the  case  of  zero  yarn/yarn  and  projectile/fabric  friction  conditions 
are  presented  and  discussed. 

Examples  of  the  temporal  evolution  of  deformation  within  the 
fabric  obtained  using  the  yarn-level  and  the  unit-cell  based  FEM 
analyses  in  the  absence  of  yarn/yarn  and  projectile/fabric  friction 
are  displayed  in  Figs.  8a-d  and  9a-d,  respectively.  In  these  figures, 
side  and  top  views  of  the  fabric  along  with  superimposed  contour 
plots  of  the  transverse  displacement  are  shown.  A  simple  compari¬ 
son  of  the  corresponding  results  reveals  that  the  temporal  evolu¬ 
tions  of  the  deformation  state  of  fabric  are  quite  similar  in  the 
two  FEM  analyses  at  shorter  post-impact  times,  Figs.  8a  and  d  and 
9a  and  d,  respectively.  At  longer  post-impact  times,  Figs.  8c  and  d 
and  9c  and  d,  the  agreement  is  less  satisfactory.  The  main  reason 
for  the  observed  discrepancy  is  that  in  the  yarn-level  FEM  analysis, 
the  principal  yarns  can  readily  move  sideways  to  avoid  the  projec¬ 
tile,  while  such  yarn  mobility  is  not  present  in  the  unit-cell  based 
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Fig.  9.  The  temporal  evolution  of  deformation  in  the  fabric  for  the  unit-cell  based  FEM  model  under  the  yarn/yarn  /ry/y  and  projectile/fabric  /rp//  friction  coefficients  of  0.0. 
Contour  bands  correspond  to  different  values  of  the  transverse  displacement,  i.e.  the  displacements  normal  to  the  fabric  surface.  All  four  fabric  edges  are  fixed. 


FEM  analysis.  Consequently,  less  yarn  breaking  takes  place  in  the 
former  case,  and  after  the  projectile  has  fully  penetrated  the  armor, 
the  “pushed  aside ”  yarns  tend  to  recoil.  The  yarn  recoiling  process 
reduces  both  the  size  of  the  penetration  hole  and  the  extent  of  resid¬ 
ual  lateral  displacements  in  the  fabric,  Figs.  8d  and  9d,  respectively. 
Despite  these  discrepancies  the  overall  final  sizes  of  the  penetration 
hole  are  comparable  in  the  two  analyses  and,  more  importantly,  the 
residual  velocities  of  the  projectile  obtained  in  the  two  analyses 
(398.6  m/s  in  the  case  of  the  yarn-level  and  397.8  m/s  in  the  case 
of  the  unit-cell  based  FEM  analysis)  are  quite  comparable. 

In  regard  to  the  effect  of  friction  on  the  fabric  deformation  and 
failure  behavior  in  the  case  of  the  yarn-level  FEM  analysis,  Figs.  6a- 
d  and  8a-d,  the  following  main  observations  are  made: 

(a)  No  significant  changes  are  introduced  initially  in  the  tempo¬ 
ral  evolution  and  the  spatial  distribution  of  the  transverse- 
deflection  wave  due  to  elimination  of  the  yarn/yarn  and 
the  projectile/fabric  friction. 


(b)  The  structure  of  the  fabric  (including  its  deformation  and 
failure)  in  the  impact  region  and  in  the  nearby  surrounding 
regions  is  greatly  affected  in  the  case  of  zero  yarn/yarn  and 
projectile/fabric  friction  conditions.  Specifically,  in  the 
zero-friction  case,  yarns  were  substantially  displaced  in  the 
in-plane  directions  away  from  the  center  of  impact.  This 
finding  can  be  readily  explained  by  the  fact  that  the  friction 
at  the  yarn  crossovers  provides  resistance  to  the  relative  tan¬ 
gential  motion  of  the  yarns,  while  such  resistance  is  absent 
in  the  zero-friction  case.  Consequently,  in  the  zero-friction 
case  yarns  impacted  by  the  projectile  are  pushed  outward, 
a  fewer  number  of  yarns  are  broken  and  the  projectile  man¬ 
ages  to  penetrate  the  fabric  mainly  by  “ wedging ”  through  it. 

A  final  examination  of  the  corresponding  results  displayed  in 
Figs.  8a-d  and  9a-d  suggests  that,  considering  the  fact  that  the 
FEM  mesh  is  significantly  finer  and  can  capture  more  details  in 
the  case  of  the  yarn-based  FEM  analyses,  the  overall  agreement  be- 
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Fig.  10.  The  effect  of  yarn/yarn  nyjy  and  projectile/fabric  fip/f  friction  coefficients  on 
the  temporal  evolution  of  projectile’s  kinetic  energy  loss  and  on  different  energy¬ 
absorbing  mechanisms  in  the  fabric  obtained  using  the  yarn-level  FEM  model:  (a) 
Hy/y  =  f^pif=  0-5  and  (b)  fiyiy  =  /^p//=  0.0. 

tween  the  two  analyses  relative  to  the  extent  of  deformation  and 
damage  in  the  armor  can  be  deemed  as  less  satisfactory  than  in 
the  case  when  friction  is  present  at  the  yarn/yarn  crossovers  and 
projectile/fabric  contact  surfaces,  Figs.  6a-d  and  7a-d.  The  results 
displayed  in  Figs.  6-9  were  obtained  under  the  projectile  velocity 
of  400  m/s.  A  similar  level  of  agreement  between  the  correspond¬ 
ing  results  obtained  in  the  two  FEM  analyses  was  obtained  when 
the  projectile  velocity  was  doubled  or  halved  relative  to  this  value. 

4A.3.  Fabric’s  energy-absorption  potential 

In  this  section,  temporal  evolutions  of  the  three  main  energy¬ 
absorbing  mechanisms  (i.e.  yarn  strain  energy,  yarn  kinetic  energy, 
and  frictional-sliding  losses)  are  compared  for  the  two  types  of 
FEM  analysis  under  the  two  types  of  frictional  conditions.  The  cor¬ 
responding  results  for  the  two  types  of  FEM  analyses  are  displayed 
in  Figs.  10a  and  b  and  11a  and  b,  respectively.  In  Figs.  lOa-d  and 
lla-d,  the  quantities  displayed  along  the  vertical  axis  are  all  nor¬ 
malized  by  the  corresponding  maximum  projectile  energy  loss. 
Also,  the  residual  projectile  velocities  are  indicated  and  their  com¬ 
parison  suggests  a  quite  good  agreement  between  the  correspond¬ 


Time, micro-sec 


Fig.  11.  The  effect  of  yarn/yarn  fiy/y  and  projectile/fabric  /^//friction  coefficients  on 
the  temporal  evolution  of  projectile’s  kinetic  energy  loss  and  on  different  energy¬ 
absorbing  mechanisms  in  the  fabric  obtained  using  the  unit-cell  based  FEM  model: 
(a)  py,y  =  iip/f=  0.5  and  (b)  fiyly  =  iiplf=  0.0. 

ing  values  obtained  using  the  two  analyses.  A  simple  comparison  of 
the  results  displayed  in  Figs.  10a  and  b  and  11a  and  b  reveals  that: 

(a)  A  somewhat  larger  contribution  of  the  frictional-sliding  loss 
is  present  in  the  case  of  the  yarn-level  based  analysis, 
Fig.  10a  vs.  Fig.  11a.  This  finding  is  consistent  with  the  fact 
that,  in  the  case  of  the  unit-cell  based  analyses,  yarn/yarn 
friction  is  not  included  explicitly  and  does  not  contribute 
to  the  frictional-sliding  loss. 

(b)  The  yarn  strain  energy  contribution  for  both  friction  condi¬ 
tions  is  higher  in  the  case  of  the  unit-cell  based  analysis. 
While  this  is  expected  considering  the  fact  that  the  effect 
of  yarn/yarn  friction  is  included  implicitly  in  the  unit-cell 
based  model  through  the  friction-dependent  in-plane  shear 
term,  the  same  argument  cannot  be  made  in  the  zero-fric¬ 
tion  case.  Thus  the  unit-cell  based  analysis  tends  to  some¬ 
what  over-predict  the  strain  energy  contribution. 

(c)  In  the  case  of  the  yarn-level  FEM  analysis,  the  full-penetra¬ 
tion  time  (defined  as  a  time  when  the  projectile  reaches  its 
final  velocity)  is  greater  (~26  ps)  in  the  presence  of  friction, 
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Fig.  10a,  than  (~15  jis)  in  the  zero-friction  case,  Fig  10b.  This 
finding  is  related  to  the  fact  that  in  the  absence  of  friction 
the  projectile  can  wedge  through  the  fabric,  while  in  the 
presence  of  friction  more  yarn  failure  is  required  before 
the  projectile  can  reach  its  residual  velocity.  In  the  case  of 
the  unit-cell  based  FEM  analysis,  the  presence/absence  of 
friction  has  less  effect  on  the  penetration  time  since  no 
“ wedge  through ”  effect  can  take  place,  Fig.  11a  and  b. 

(d)  The  contribution  of  the  frictional-sliding  mechanism  to 
energy-absorption  is  typically  small.  However,  the  presence 
of  friction  itself  promotes  projectile’s  kinetic  energy-absorp¬ 
tion  and,  thus,  increases  fabric  energy-absorption  capacity 
by  yarn  straining  and  yarn  acceleration.  This  finding  is  based 
on  a  comparison  between  the  residual  projectile  velocities  in 
the  friction  case,  Figs.  10a  and  11a  and  in  the  zero-friction 
case,  Figs.  10b  and  lib. 

(e)  As  the  projectile  begins  to  interact  with  the  fabric,  the  con¬ 
tributions  of  yarn  straining  to  the  energy-absorption  process 
increases.  Once  the  projectile  starts  to  penetrate  the  armor, 
further  yarn  straining  either  ceases  or  increases  at  a  signifi¬ 
cantly  lower  rate.  After  the  armor  is  fully  penetrated,  the 
fabric  kinetic  energy  component  of  the  total  yarn  energy 
experiences  a  relative  increase  (due  to  yarn  recoiling)  and 
likewise  the  frictional-sliding  component  of  the  total  yarn 
energy  increases. 

(f)  The  extent  of  agreement  between  the  two  FEM  analyses 
regarding  the  energy-absorption  capacity  of  armor  under 
the  projectile  velocity  of  400  m/s  as  seen  in  Figs.  10a  and  b 
and  11a  and  b  was  also  observed  under  different  initial  pro¬ 
jectile  velocities. 


4.2.  The  effect  of  far-field  boundary  conditions 

The  results  presented  and  discussed  in  Section  4.1  pertain  to  the 
case  of  fixed  boundary  conditions  applied  to  all  four  edges  of  the 
armor.  It  is  generally  found  that  the  nature  (free  vs.  fixed)  of 
boundary  conditions  affects  the  ballistic  response  of  the  armor 
and  its  ability  to  absorb  kinetic  energy  as  well  as  it  affects  the  rel¬ 
ative  contributions  of  various  energy-absorbing  mechanisms  (e.g. 
[10-13]).  The  results  presented  and  discussed  in  Section  4.1  show 
that  the  meso-scale  unit-cell  based  material  model  can  reasonably 
well  account  for  the  deformation/failure  behavior  and  the  energy¬ 
absorbing  capacity  of  a  plain-woven  single-ply  fabric  under  fixed 
boundary  conditions  being  applied  to  all  four  fabric  edges.  In  this 
section,  a  comparison  is  made  between  the  corresponding  results 
obtained  using  the  yarn-level  FEM  and  unit-cell  based  FEM  analy¬ 
ses  under  other  boundary  conditions.  Due  to  space  limitations,  the 
results  analogous  to  those  displayed  in  Figs.  6-11,  but  for  the 
boundary  conditions  corresponding  to:  (a)  two  free  (opposite) 
and  two  clamped  (opposite)  fabric  edges,  and  (b)  four-free  fabric 
edges  will  not  be  shown  in  this  section.  Instead  these  results  will 
be  discussed  and  related  to  their  counter  parts  presented  in  the 
previous  section. 

Regarding  the  spatial  distribution  and  temporal  evolution  of  the 
transverse  displacements  within  the  fabric,  similar  observations 
were  made  as  those  related  to  Figs.  6-9.  The  main  differences  were 
observed  near  the  edges  of  the  fabric  along  which  the  boundary 
conditions  were  changed  from  clamped  to  free.  Specifically,  when 
the  transverse  wave  front  reaches  the  free  fabric  edges,  the  edges 
begin  to  be  pulled-in  towards  the  point  of  impact  resulting  in 
“ bowed-in ”  edges.  This  situation  is  similar  to  the  one  observed  in 
sheet-metal  forming  where  reduction  in  the  clamping  force  may 
lead  to  a  transition  from  stretch-forming  to  deep-drawing.  As  a  re¬ 
sult  of  this  transition,  in  the  present  case  the  largest  maximum 


transverse  deflections  were  obtained  under  the  four-free  edge 
boundary  conditions. 

The  agreement  between  the  extents  of  fabric  damage  predicted 
by  the  two  FEM  analyses  under  the  two  types  of  boundary  condi¬ 
tions  analyzed  in  this  section  is  at  a  level  similar  to  that  observed 
in  the  previous  section. 

Regarding  the  energy-absorbing  capacity  of  the  fabric,  it  is  ob¬ 
served  that  this  capacity  is  increased  when  the  boundary  condi¬ 
tions  applied  along  the  edges  of  the  fabric  are  changed  from 
clamped  to  free.  This  observation  is  consistent  with  the  fact  that 
free  boundary  conditions  lower  the  rate  of  stress  increase  in  the 
fabric  and,  thus,  delays  fabric  failure.  Consequently,  the  fabric  is  gi¬ 
ven  longer  time  to  interact  with  the  projectile  and  absorb  its  ki¬ 
netic  energy  before  penetration. 


5.  Summary  and  conclusions 

Based  on  the  material  model  development  procedure  utilized 
and  the  results  of  the  subsequent  computational  analyses,  the  fol¬ 
lowing  main  summary  remarks  and  conclusions  can  be  drawn: 

1.  A  simple  meso-scale  unit-cell  based  material  model  for  plain- 
woven  single-ply  fabric  armor  is  developed  and  implemented 
in  a  User-material  Subroutine  suitable  for  use  with  commercial 
explicit  finite  element  programs. 

2.  The  User-material  Subroutine  is  coupled  with  the  ABAQUS / 
Explicit  finite  element  program  and  used  in  a  series  of  transient 
non-linear  dynamic  analyses  of  a  plain-woven  single-ply  fabric 
armor  impact  by  a  spherical  steel  projectile  in  order  to  test  and 
validate  the  material  model. 

3.  By  comparing  the  results  obtained  in  a  finite  element  analysis  in 
which  the  User-material  Subroutine  is  used  with  those  obtained 
in  a  computationally  more-costly  FEM  model  in  which  yarn 
weaving  is  accounted  for  explicitly,  it  is  found  that  the  present 
meso-scale  unit-cell  material  model  for  plain-woven  single-ply 
fabric,  under  fixed  boundary  conditions  applied  to  all  four  edges 
of  the  fabric,  can  reasonably  well  account  for  the  observed  tem¬ 
poral  evolution  and  spatial  distribution  of  transverse  waves  and 
damage  within  the  armor  and  for  the  roles  of  various  energy¬ 
absorbing  mechanisms.  Similar  observations  were  made  where 
different  combinations  of  clamped  and  free  boundary  condi¬ 
tions  were  applied  along  the  edges  of  the  fabric. 
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